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GEOS-C ALTIMETER ATTITUDE BIAS ERROR CORRECTION 


INTRODUCTION 

A pulse-limited split-gate-tracking radar altimeter has been flown on 
Skylab and will be used aboard GEOS-C [1] . If such an altimeter were to em- 
ploy a hypothetical isotropic antenna, the altimeter output would, of course, be 
independent of spacecraft orientation. To reduce power requirements, however, 
the gain of the altimeter antenna proposed has been increased to the point where 
its beamwidth is only a few degrees. The gain of the antenna consequently varies 
somewhat over the pulse- limited illuminated region of the ocean below the al- 
timeter, and the altimeter output varies with antenna orientation. The error 
thereby introduced into the altimeter data has been modeled empirically [ 2] , 
but because of an unfortunate choice of the form of the modelling function, close 
agreements with the expected errors was not realized. In this report the attitude- 
error effects expected with the GEOS-C altimeter are more closely modelled 
using a form suggested by an analytical derivation. As in reference [ 2] , the 
treatment here is restricted to the case of a relatively smooth sea, where the 
height of the ocean waves are small relative to the spatial length (pulse duration 
times speed of light) of the transmitted pulse. 


Expected Return Power 

The expected value of the received power has been given by Harger [3] . 
Referring to Figure 1, this power is 
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Expected value of received power at time t. 

Wavelength = 2.158 cm for GEOS-C altimeter 
Radar cross section of ocean surface per unit area 
Satellite height above mean sea level = 843 km for GEOS-C 
Power of transmitted pulse as a function of time 
2r 0 /c = Roundtrip delay along the vertical* 

1. + ro /r = Earth curvature correction 
Earth radius = 6378 1cm 
Speed of light 

Antenna gain (assumed symmetric about boresight) as a function 
of the angle 'P off boresight 

Probability density function of the departure h of the sea surface 
from mean sea level. 


In deriving (1) it was assumed that the scattering is spatially incoherent 
and that the scattering region from which the received power is returned is 
confined to values of 6 of less than a few degrees. Over this small angular 
region cr 0 is taken to be constant. With the additional restriction that the sea 
wave heights are small relative to the pulse length, the first integration in (1) 
can be performed to give 
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There is no loss in generality in taking the line of boresight of the antenna to 
lie in £he</> = 0 plane at an angle with the vertical. The angle 0 between 
boresight and the line r is then given by 


cos \p = cos $ cos 9 e + sin 0 sin cos 
Neglecting atmospheric effects on propagation, which can be corrected for seperately. 
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Because of the small values assumed for $ and 8 e this becomes 


i/r 2 = d 2 + - 266 e cos 4> 


( 4 ) 


The expected value of the power returned from an incoherently-scattering 
smooth curved ocean surface may be calculated from equations (2) and (4) if the 
shape P T of the transmitted pulse and the gain function G of the antenna are 
known, the calculation being accurate provided that the illuminated area involved 
in the integration does not extend over more than a few degrees from the vertical. 


GEOS-C Altimeter Split Gate Tracking 

The GEOS-C altimeter will employ two gates for tracking [1]« These gates 
are each one pulse- width wide, and are separated by one pulse- width. The 
position of the gates is automatically adjusted so that the power in the first is 
half that in the second. This position then provides a measure of altimeter 
height, with the time difference between the start of the transmitted pulse and 
the start of the first gate taken to be the round trip delay time along the vertical. 
This operation is described by the equation 
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Where r is the roundtrip delay time measured by the altimeter and r is the 
pulse width. If the bias error in this measurement is denoted by r e , then 
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T 0 + T e 


(6) 


From (2), (5) and (6) 
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To solve equation (7) above for the bias error r e , the pulse shape P T and the 
antenna gain function G(\p) must be known. Here we assume a rectangular pulse 
of length r 


P T (x) = P Q ; 0 < x < r 


= 0; otherwise 


( 8 ) 


and, for the antenna, a circular aperture with an illumination that results in 
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where 


a 


7TB 

~ 


s in 0 


a « diameter of aperture = 60.96 cm 

Equation (7) was solved for the time delay bias r e and the corresponding 
error A r 0 = c r # /2 in the altimeter height measurement. The left hand side 
of (7) is a function of r e whose derivative can be found by means of the usual 
rule for the differentiation of an integral whose limits depend on a parameter. 
The solution was consequently readily obtained by numerical integration and 
the use of the Newton-Raphson method [4]. The APL [5] programs used and 
a sample calculation using a 200 ns pulse and a 0.1 degree pointing error are 
given in the Appendix. The bias errors Ar 0 for pulse widths of 12.5 and 200 
nanoseconds are plotted as solid lines in Figures 2 and 3 respectively. 


Fitting of Error Curves 

In order to use the curves of Figures 2 and 3 in the processing of GEOS-C 
data, it is desirable to fit the curves with a suitable function. Such a function 
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could be found by trial and error, but here we have preferred to be guided by 
an approximate analytic solution to equation (7). First the square of the gain 
(9), ■which is an even function of y> , was approximated by 


G 2 0/0 = G 2 0 
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where 


=i ^ = 1.22° = 0.0213 radians 
2 a 

is the half beam width. The constants A and B obtained by least squares were 
about 1.25 and 0.5 respectively. Substituting (10) into (7), the integration may 
be performed directly. The resulting expression is a quartic equation in r e . 
For small values of r,(r e /r<< 1), the bias error is 
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An inspection of (11) leads to the conclusion that the altitude error Ar 0 i s quite 
sensitive both to the pulse width and, through the dependence on the values of 
the constants A and B, to the exact shape of the antenna beam. An examination 
of the magnitude of the terms in (11) suggests the use of the form 


Ar 0 = - aj 





( 12 ) 
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to fit the solid curves of Figures 2 and 3. Using the values a 2 = 1.066 and a 2 =1.348 
obtained by least-squares fitting, the errors given by (12) are plotted as dashed 
lines on the Figures. Also shown on the Figures are curves calculated from the 
equation of reference 2. 


Summary 

The GEOS-C altimeter bias caused by antenna pointing error has been 
calculated assuming a calm ocean, a rectangular pulse shape and a specific 
antenna gain. In the long pulse mode, the bias error is significant, varying over 
a range of 4 meters as the pointing error varies from 0 to 1.1 degrees. A 
simple formula is given for correcting this error when the spacecraft attitude 
is known. 

This analysis represents an ongoing NASA-GSFC effort in support of the 
GEOS-C altimeter capability which has been developed under the direction of 
the NASA -Wallops Flight Center. Mr. H. Ray Stanley of the Applied Sciences 
Directorate is the GEOS-C Project Scientist. Mr. C. L. Purdy of the Engineer- 
ing Directorate, Wallops Flight Center is responsible for GEOS-C altimeter 
hardware development, data collection, data preprocessing, calibration and data 
distribution. References 6 and 7 document previous analyses related to the 
work in this report. This ongoing study is being extended to address sea state 
effects on satellite altimetry and the results will be published in a subsequent 
report. 
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Appendix. APL Programs 
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Figure 3. Altitude Error vs. Attitude Error, r = 200 ns 
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